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ABSTRACT 

The effect of magnetic fields on the frequencies of toroidal oscillations of neu- 
tron stars is derived to lowest order. Interpreting the fine structure in the QPO 
power spectrum of magnetars following giant flares reported by Strohmayer and 
Watts (2006) to be " Zeeman splitting" of degenerate toroidal modes, we estimate 
a crustal magnetic field of order 10 15 Gauss or more. We suggest that residual m, 
-m symmetry following such splitting might allow beating of individual frequency 
components that is slow enough to be observed. 

Subject headings: stars: magnetic — pulsar: neutron — stars: oscillations — X-rays: 
stars 



The disco very of quas i -periodic oscillations (QPO) in th e hyperflares of soft gamma- 



ray repeaters (jlsrael et al.l l2005t IStrohmayer fc Watts 1 120051 ; IWatts fc Strohmayer 1 12006 



Strohmayer fc Watts 1 120061 ) has attracted much attention to the study of nonradial os- 
cillations of neutron stars wit h solid crust. These oscillations w ere studied extensively 



in ( 


Hansen & Cioffi 


1980: 


Unno et al. 


( 


19891) discuss 



Jl989h discuss in their book many aspects of the theory of nonradial oscilla- 



tions of stars. The properties of non radial modes of strongly magnetized neutron stars have 



been invest igated by several authors (lDuncanlll998l ; iPiro 1120051 ; iLee 1120071 . 120081 ; ISotani et al. 



20061 . 120071 ). Their main focus was on the study of axisymmetric modes. Duncan anticipated 
toroidal oscillations as a result of giant flares and noted that the magnetic field could affect 
their frequencies. The effects of a strong ver tical magnetic field on the oscillation spectrum 
of a cylindrical slab model were studied in (jCarroll et al.lll986l ). The nonradial modes are 
generally divided into two main classes: the spheroidal and toroidal modes. We will concen- 
trate on the study of toroidal modes because of their possible connection with QPOs. They 
are defined by conditions 
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V -u = 



where u is a displacement vector. 



Without magnetic field and rotation, toroidal modes ( denoted it n , where the index n 
is the number of radial nodes in the eigenfunction ) have frequencies that do not depend on 
m. This degeneracy is lifted by magnetic field since it breaks the spherical symmetry of the 
problem. In this paper we will study how this happens considering the magnetic field B as a 
perturbation. T he influence of a m agnetic perturbati on on spheroidal modes was considered 
earlier [see e.g. (jUnno et al.l 119891 )]. We also refer to iDahlen fc Trompl (119981 ) for a detailed 
discussion of perturbation theory with applications in seismology. Assuming an oscillatory 
time dependence u oc e~ lujt , where u is the mode frequency, the equation for eigenfunctions 
without a magnetic field is 

- uj 2 u = A(u) (2) 

where the linear operator A describes the dynamics given specific parameters for the neutron 
star. We do not need its exact form here. Inte rested readers may consult e.g. ([Hansen fc Cioffi 
19801 ; iMcDermott. van Horn &: Hansenlll988l ). Let us denote toroidal eigenfunctions and fre- 
quencies obtained from this equation as ul; and In spherical coordinates components 
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Here w 



(o) 
i 



is a radial eigenfunction. 



With a magnetic field we have 



uj 2 u = A(u) + ^r[(V x b) x B) 



where 



b = V x (ux B] 
We now take B to be the uniform field 

B = Bne*. 



(4) 
(5) 

(6) 



It is helpful to use vector spherical harmonics [see e.g. (jVarshalovich et al.lll988l )] to express 
the vector operators in equations (03) and ([5]) in Y\ m representation. Definitions of vector 
spherical harmonics and some usef ul formulae a re giy en in the appendix. Then, using the 
perturbation approach described in lUnno et al.l (Il989l ). we obtain 



Ulr, 



i ; (0) + , 



(7) 
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where 



LO 



(i) 

Im 



(o) 



LO 



where 



(0)2 r (UJ2 o j 



(0) 




F(l,m) 



/(/ + 2)((/ + l) 2 -m 2 ) 



m 



+ 



(Z 2 - 1) (Z 2 - m 2 ) 



F(Z,m) (8) 



(9) 



(Z + l) 2 (2Z + l)(2Z + 3) Z 2 (Z + 1) 2 ' Z 2 (2Z + 1) (2Z - 1) 
Table 1 lists values of F(l,m) for specific modes. 

This is the main result of our paper within the framework of a perturbative calculation. 
Perturbation theory is applicable when 



LO 



(i) 
i in 



LO 



(0) 



< 1 



(10) 



For the case of a spherical star with uniform shear modulus /i, one can calculate the integral 
in (jSJ) exactly and obtain 



(i) 

im 



LO 



(0) 



87T/1 



F(l,m) 



(11) 



By ffTD]) perturbation theory is thus applicable when 



3 

87T/1 



< 1 



(12) 



Because of the axisymmetry of the magnetic field the (21+1) degeneracy in m is only partially 
lifted : it depends only on \m\. The mode splits into (1+1) modes. 

The "92 hertz" mode in the QPO data following the 2004 giant flare of SGR 1806-20 
actually seems to be accompani ed by a conglomerate of man y peaks in the power spectrum , 
ranging from about 78 to 105 hz [Strohmayer fc Watts I (120061 ). IWatts fc Strohmayer 1(120061 )]. 

implying that at it most extreme value -fey ~ 0.25. The 92 hertz QPO is attributed to 



the 1=7 toroidal mode jStrohmayer &: Watts I (120061 ). IWatts fc Strohmayer I (120061 )]. There 
also seems to be a significan t component at about 80 hz [panels 1-4, and 17 of figure 9 in 
Strohmayer & Watts I (120061 )]. which could be the 1=6 mode. [If the field is axisymmetric 
and has mir ror symmetry around the equator , then the mechanism for luminosity variation 
proposed by iTimokhin. Eichler &: Lyubarsky I (120071 ) works for odd 1 modes, where the two 
foot points of a given magnetic field line move in opposite directions. This introduces a twist 
in the magnetic field line that implies a current perturbation. In the more likely situation 
that the field lacks this high degree of symmetry, even 1 modes are also possible.] However, 
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there are also modes in the 80-90 hz range and also significantly above 92, up to 100 hz 
or more, so it is hard to interpret the data unambiguously. While there is some evidence 
for systematic increase of the frequency with time, the QPO signal in any given time frame 
appears to be constant, and in some cases several bands appear simultaneously, separated 
by several hz. We conclude that if the splitting is due to a magnetic field, it is at least 2 
percent, and at most about 25 percent. 

B 2 

According to the calculations here, the variation is expected to be of order 0.4 to 0.5g^, 

as m ranges from to 7. This implies that ~ 0.04, or B ~ 0.2a/87t/i if we interpret the 
2 percent splitting to be due to magnetic effects. 

If the entire range of from 78 to 105 hz is attributed to magnetic splitting, our first 
order approximation is only marginal at this strength, and a higher order calculation would 
give a slightly lower value for the frequency shift by a given field strength, so all we can say 
is that the field is of the same order as (Sufi) 1 / 2 . 



The value of (Aufi) 1 ^ 2 has been estimated by to be about 6x 10 15 Gauss ( Thompson. C. &: Duncan. R. C 



( 119951 )). so, if we assume a magnetic splitting of 2 percent, the value of B appears to be of 
order 1.7 x 10 15 Gauss or higher, i n reasonably good a g reement with that es timated for the 



dipole field component, 1.6 x 10 15 (IWoods et all (120021 ): IPalmer et all (120051 )1 Note that in 



the geometry used here, the field lines do not particularly lie within the crust but rather cut 
through it vertically at angle 7r/2 — 9, where 9 is the latitude, so a purely toroidal field would 
be estimated to be somewhat weaker than the above estimate by a factor of 30 percent or 
so. On the other hand, a toroidal field could easily be somewhat larger than the poloidal 
field without affecting the dipole moment. 

One possible observational consequence of magnetic splitting of frequency degeneracy 
for toroidal oscillations is the fact that it leaves m and -m modes degenerate to the extent 
that the field is axisymmetric about the magnetic axis. In contrast to the magnetic splitting 
between different \m\ modes (which appears to be on the order of several hertz and implies 
that the beat periods would be less than the rotation period), the m and -m modes would 
beat more slowly, and their beating could possibly be observed on a timescale of perhaps 
seconds to tens of seconds. The appearance of a particular frequency band beating on 
and off during the QPO activity could be a signature of the m and -m modes having [or, 
more precisely, of their symmetric and antisymmetric combinations having] slightly different 
frequencies. Such beating would be a measure of a non-axisymmetric component to the field. 
If, for example, the field is larger at <fi = than at <fi = tt/2, then the symmetric combination 
of the m=l and m=-l modes, proportional to cos0, would have a slightly higher frequency 
than the anti-symmetric mode, proportional to sin0. 
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The residual m,-m degeneracy could also be removed by rotation of a neutron star. For 
toroidal modes, the angular frequency shift d ue to rotation i n a r ot ating referen c e fram e 



attached to the star is mVt/l{l + 1) (see eg. iPekeris et al. I (119611 ); IStrohmayer I (119911 ) ) . 
where Q is the angular frequency of rotation. With period of rotation P rot ~ 7.5s one 
can see that the splitting is quite small and modulations of the crustal displacement am- 
plitude (beats) are possible with period T mo d = P ro tl(l + I) /2m. Since the physics of 
QPO variations in magnetar lum i nosity is somehow determined by crustal oscillations [e.g. 



Timokhin. Eichler fc Lyubarsky I (120071 )]. one might also expect time modulations in the 



observed QPO components on time scales of order 30 seconds at I = m = 7. 

Another effect that could be looked for is the different damping rates of the symmetric 
and antisymmetric combinations due to coupling with the core. In general, such damping 
occurs when the frequency of the mode matches the resonant Alfven frequency of some 
connecting field lines in the core. But the damping rate would depend on the amplitude of 
the crustal oscillation at the latitude and longitudes where the frequency match happens to 
take place, and this amplitude can differ among the various linear combinations of symmetric 
modes. This suggests that QPO components could resolve to narrower frequency bands as the 
more rapidly damped linear combination gives way to the surviving combination. Detecting 
such an effect, however, would require good frequency resolution. 

The differences between beating, resolution to the longest lived of several modes, and 
continuous frequency drift should all be made clear. We are unable to see how magnetic 
splitting or damping by the continu um leads to continuou s frequ ency drift, and we see little 



if any evidence for it in the data of IStrohmayer fc Watts I (120061 ). 



The potential wealth of data available in magnetar seismology awaits confirmation of 
an accepted model for it. Future observations of intermediate flares, which may be more 
frequent than giant flares, may provide badly needed additional data. 

The authors thank Drs. A. Watts and Y. Lyubarsky for helpful discussions. We ac- 
knowledge the Israel-U.S. Binational Science Foundation, the Israel Science Foundation, and 
the Joan and Robert Arnow chair of Theoretical Astrophysics for generous support. 



A. Vector spherical harmonics 

Vector spherical harmonics are defined as 

YjM C t^ Y Lm (B, <p) e a 
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where 



1 ^ • —> \ — * - - 1 ^ 
J+i = -—^ [e x + «e y ), e = e z , e_ : = -j= (e a 



Z 6V 



Here L can have values L = J, J ± 1 for a given J . 

During the calculati ons following formulae prove useful (for more information see e. 
JVarshalovich et al.l[l988h \ 
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Now, taking into account that and e 2 can be represented as [e.g. f Varshalovich et al 
19881 )] 
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one can use the above formulae to get the result (IH1I91) . 
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Table I. Values of the function F(l, m) for I — 5, 7, 9 



m 


F(5,m) 


F(7, m.) 


F(9,m) 





0.487 


0.493 


0.496 


1 


0.472 


0.485 


0.49 


2 


0.426 


0.459 


0.474 


3 


0.349 


0.416 


0.447 


4 


0.241 


0.356 


0.409 


5 


0.103 


0.279 


0.36 


6 




0.184 


0.3 


7 




0.074 


0.23 


8 






0.149 


9 






0.057 



